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ABSTRACT. For a single-server multi-station polling system we focus on the generating function and 
Laplace-Stieltjes transform of the time-stationary joint queue length and workload distributions, re- 
spectively, under no further assumptions on the service discipline. We express these quantities as 
expressions involving the generating functions of the joint queue length distribution at visit begin- 
nings and visit completions at the various stations. The latter is known for a broad variety of cases. 
Finally, we identify a workload decomposition result. 



1. Introduction 

Consider a queueing model consisting of multiple queues: Qi, . . . ,Qn, N > 1, attended by a sin- 
gle server, visiting the queues one at a time in a cyclic order. Moving from one queue to another, 
the server incurs a (possibly zero) switchover time. Such single-server multiple-queue models are 
commonly referred to as polling systems. Triggered by a wide variety of applications, polling sys- 
tems have been extensively studied in the literature; see [11, 14, 16] for a series of comprehensive 
surveys and [1, 8, 12] for extensive overviews of the applicability of polling systems. 
In the polling literature much attention has been given to the determination of the probability 
generating function (PGF) of the joint queue length distribution of the polling system at various 
epochs. In particular, the PGF's of the joint queue length distribution at the epochs at which the 
server begins a visit, % ( (■), and completes a visit, ^ (•), at Qi were extensively studied. A seminal 
result has been obtained by Resing [9], who found these PGF's for polling systems with service 
disciplines satisfying the so-called branching property. Service disciplines - rules which establish 
the server's behaviour while visiting a queue - that satisfy the branching property include the 
well known exhaustive discipline (per visit the server continues to serve all customers at a station 
until it empties) and gated discipline (per visit the server serves only those customers at a station 
which are found there upon its visit). 

Polling systems with disciplines which do not satisfy the branching property usually defy exact 
analysis. A typical example of a discipline that does not satisfy the branching property is the k- 
limited discipline (per visit at most k customers are served at a station). Nevertheless, it is possible 
to compute % 4 and Y Ci for some special cases like completely symmetrical queues or like the two- 
queue case. For example, for the case N = 2 with Qi served exhaustively and Q2 according to the 
fc-limited discipline, the formulas for % i and ~f Ci were found by Winands et al. [17]. Furthermore, 
for the case of TV = 2 and Q\, Q2 both served according to the 1-limited discipline, the formulas 
for and "¥ Ci were found by solving a boundary value problem [4]. For other examples of 
polling systems with non-branching disciplines that admit a solution for these PGF's we refer to 
the surveys [13, 14]. 
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Although the PGF's % i and y Ci have been found in a very broad variety of cases, the PGF 
of the steady-state joint queue length distribution at an arbitrary epoch, which in principle is the 
most important PGF, received little attention and has not been studied in full generality, one ex- 
ception (that we are aware of) being Sidi et al. [10]. The first contribution of the present paper is 
to determine M for a quite general class of polling systems. To explain the formula for =S let us 
introduce two additional PGF's. Let 5?^ (•) and (•) be the PGF's of the joint queue length dis- 
tribution at service beginnings and service completions, respectively. Eisenberg [7] established a 
link between the various PGF's: % u "f Ci , S*b- and 5? Ci . In fact, from his observation it follows that 
it is enough to determine, for example, i = 1, . . . , N, to know all the remaining PGF's y ci , S% it 
J^ Ci , i = 1, . . . , N. Hence, as noted before, all these PGF's are well known for a vast class of polling 
systems. One of our main results - Theorem 1 - reveals that the formula for i? can be expressed 
in terms of ,5^ Ci , i = 1, . . . , N (or equivalently in any other aforementioned PGF's). Using a similar 
technique that allowed us to determine we were also able to find the Laplace-Stieltjes 

transform (LST) of the joint workload distribution in steady state. We give a formula for W in 
Theorem 2 expressed in % i and f Ci , which is the second contribution of our paper. 
J? was found in Sidi et al. [ 1 ] for polling systems restricted to the case of cycles with either 
all disciplines being gated or all being exhaustive. W was found in [ ], for the class of polling 
systems with Levy input, but restricted to the case of branching-type service disciplines. Both 
of the results from [5, 10] assume non-zero switchover times. It should be emphasized therefore 
that our main results - Theorem 1 and Theorem 2 - hold regardless of the service discipline and 
provide a unified expression for £2 and W for both the zero and non-zero switchover case; see 
also Remark 1 and Remark 4. 

From Theorem 2 we were also able to retrieve the well known decomposition of the total workload 
W in a polling system into an independent sum of the amount of work Wm/g /i m a corresponding 
M/G/l queue and W SW it c h, the total amount of work in the polling system at an arbitrary epoch 
in a switching period. This decomposition was found in [3], however the distribution of W swltc h 
was only known in a few cases, cf. [15]. In Remark 5 we present the LST of the distribution of 
W r SW itch expressed in % i and Y Ci . This constitutes the third and final contribution of our paper. 
The paper is organised as follows. In Section 2 we give a detailed description of the polling model 
under consideration, introduce the various PGF's ^ , ^ , S% , and relate them to each other. 
Section 3 contains our first main result, Theorem 1 - the expression for £}, the PGF of the joint 
queue length distribution in steady state at an arbitrary epoch, in terms of -5f Ci . Section 4 contains 
our second main result, Theorem 2 - the expression for W , the LST of the joint workload distri- 
bution in steady state at an arbitrary epoch, in terms of 7^ and "V Ci . This section also contains 
Remark 5, which establishes the workload decomposition in polling systems. 



2. Model description 

2.1. A general polling system. We consider a system of N > 1 infinite-buffer queues Qi, ■ ■ ■ ,Qn 
and a single server S. The service times of customers in Qi are i.i.d. (independent, identically 
distributed) positive random variables generically denoted by B t . We denote the LST of Bi by 
Bi(-) and assume that bi := KBi < oo. The server moves among the queues in a cyclic order. When 
S moves from Qi to Qj+i, it incurs a switchover period. The durations of successive switchover 
times are i.i.d. non-negative random variables, which we generically denote by Si. We denote 
the LST of Si by Si(-) and assume that := ESi < oo; let s := Yli=i s »- Customers arrive at Qi 
according to a Poisson process with rate A;; let A := J^ili ^i- W e do not assume anything about 
the service disciplines at Qi. Define p; := Xibi as the traffic intensity at Qi', let p := Pi- m what 

follows we shall write z for an TV-dimensional vector in Mr, z = (zi, . . . , zjf), and we assume that 
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\zi\ < 1 for every i = 1, . . . ,N. Throughout the paper we implicitly use the convention that any 
index summation is modulo N, for example Qn+i = Qi- 

We assume that all the usual independence assumptions hold between the service times, the 
switchover times and the interarrival times. We assume that the ergodicity conditions are ful- 
filled and we restrict ourselves to results for the stationary situation. 

We also consider the variant of this model in which all the switchover times are zero. In the latter 
model, £ makes a full cycle (viz., passes all the queues once) when the system becomes empty 
and subsequently stops right before Q\. All this requires zero time. When the first new customer 
arrives, 5" cycles along the queues to that customer. Note that the behavior of the server when 
the system is empty does not affect the queue length distributions. However, it may involve a 
number of (zero-length) visits, and hence it does affect the queue length distribution embedded at 
visit beginning and visit completion instants. The convention that the server stops at some fixed 
queue was used in [ ] and turned out to result in cleaner expressions for that embedded queue 
length distribution than when the server is assumed to stop right away. Making a full cycle before 
stopping may be necessary if there is no central queue length information available to the server 
so that the only way to detect an empty system is to keep track if some fixed queue is visited twice 
in zero time. 



2.2. Distributional identities. For details of this subsection we refer to [ .]. 

Note that each time a visit beginning or a service completion occurs, this coincides with either 
a service beginning or a visit completion. All service beginning epochs in a visit to Qi are also 
service completion epochs at Qi, except the first service beginning epoch, which is also a visit be- 
ginning epoch. Similarly, all service completion epochs at Q t are also service beginning epochs at 
that queue, except the last service completion epoch, which is also a visit completion epoch. This 
has been observed by many authors, but most notably by Eisenberg [ ], who restricted himself 
to the cases of either exhaustive or gated service at all queues. However, its applicability is not 
restricted to a particular service discipline. After some manipulations this observation yields: 

(i) 7i y bi (z) + y Ci (z) = y bi (z) + 7i r Ci (z). 

Here % t (z) and (z) denote the PGF's of the joint queue length distribution at visit beginnings 
and visit completions at Qi, while (z) and 5^ Ci (z) denote the PGF's of the joint queue length 
distribution at service beginnings and service completions at Qi, the coefficient ji represents the 
reciprocal of the mean number of customers served at Qi per visit, i.e., the ratio of visit beginnings 
to service beginnings. With C denoting the steady-state cycle length, 1/7, = A^EC. In the case of 
non-zero switchover times, we simply have EC = s/(l — p). In the zero switchover times case 

Indeed, (0)/EC denotes the mean number of times per time unit that S arrives at Q\ to find 
the whole system empty (^ (0) is the probability that S sees an empty system at a Qi visit). Each 
of those epochs is followed by exactly one arrival of a customer to an empty system. By PASTA, 
there are on average A(l — p) such arrivals per time unit. 
We rewrite (1) into: 

(3) li {% i {z)-% i {z))=y hi {z)-y Ci {z). 

Let £(z) := Y^j=\ ^iO- ~~ z j)- ^ is easy to see that ,5^ Ci (z) and ^^(z) are related via 

(4) ^ Ci (z)=^ i( z)^M. 
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It follows from (1) and (4) that 

(5) -^W = — ¥kr^ ~ ^W) ■ 

Zi - Bi{Y,(z)) 

Next we relate to ~f Ci . In a polling system with switchover times 

(6) % i+1 (z) = % i (z)S i (Z(z)), i=l,...,N. 
In a polling system without switchover times 

(7) r bi+1 (*) = i = 1,...,JV-1. 

The relation between % ± and ~f CN deserves special attention in the zero switchover case, because 
of our convention stated in Section 2.1, concerning the behavior of the server when the system is 
empty. When all queues in the model without switchover times become empty, S in our conven- 
tion makes a full cycle and subsequently stops right before Q\ (all this requires no time). When 
the first new customer arrives, S cycles along the queues to that customer. The consequence of 
this is that when the system is empty at the start of a visit to Q\, then the next visit to Qi does not 
take place until a customer has arrived. Therefore, 



N 



(8) r bl {z) = r CN {z) - %m ( i - t z = r ^ (*) " ^t^s(«). 



i=i 



A / Ny ' A 



In particular, observe that % 1 (0) = "^c N (0)/2, due to the specific server movement in an empty 
system. 

3. Queue length distribution at an arbitrary time 

Theorem 1. For a general polling system as introduced in Section 2, let J3(-)be the probability generating 
function of the joint queue length distribution at an arbitrary time in steady-state. Then, 

(9) m , _ J_ v ( * (i - £(£(*))) i%{z)-n i+1 (z) 

W ^- EC Z^ E(z) Zi -Bi(Y,(z)) + E(*) 

ivhereEC = s/(l — p)in the non-zero switchover case and EC = % t (0)/(A(l — p)) in the zero switchover 
case. Equivalently, 

(10) * {Z) ^W-«W*\ 

E»=i M 1 - z i) 

Proof We will divide the proof into three parts. First we prove (9) for i? in the non-zero switchover 
times case. Next we will show that this formula is also valid in the zero switchover times case. 
Finally, we will show that the unified formula can be expressed as (10). 

Let Xi (-) and Yi (•) be the PGF of the joint queue length distribution at an arbitrary moment during 
a visit to Qi and during a switchover (idle) period between Qi and Qi+i, respectively. 
Non-zero switchover times: By the stochastic mean value theorem 

1 N fh 

(11) £>(z) = —J2[-Mz) + ^Y l (z) 
Furthermore, for any i, 



(12) Xi(z) = y bi (z)Br\nz)) = 3%(z) 



l-B^jz)) = ll z l {%Xz)-% l {z))\-B l {Y,{z)) 
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where i?f ast (-) is the LST of B\ asi , the past part of B 4 and where we used (5). Analogously, 



(is) Y iiz ) = %i{z) sr^z)) = y Ci ( z ) nz)Si Si nz) , 

where S* ast (-) is the LST of S\ ast , the past part of Si and where we used (6). Combining (11)-(13) 
leads to (9). 

Zero switchover times: Observe that in this case we can write 



1-Si(£(z)) _ 1 Y Cz (z) - % z+1 (z) 



B(z) = pB sermng {z) + (1 - p)B 



n on — serving 



where =S S 



,(■) and M. 



non — serving 



(•) denote the PGF of the joint queue length distribution at an 



arbitrary moment when S is serving and non-serving, respectively. Trivially, =S„on- serving ( z ) = 1- 
The stochastic mean value theorem gives 

JV 



Note that 



Hence, in the zero switchover case, 



A ? 



r hi (z)-r Ci (z) 



(l-fl,(E(z)) 



EC ^ \ E(z) 
Noting that (2), (7) and (8) give 

JV 



+ l-p. 



EC 



4=1 



£(z) 



ECS(z) 



AEC 



we conclude that (9) is also valid in the zero switchover times case. 
Simplification of (9): Observe that 

^ ^(z)-% l+1 (z) = A % % (z) - % t (z) 



A' 

E 



£(*) £r s(z) 

fir 



»=1 v ' i=l 

Hence, using (4), (5) and the definition of 7.; we have 



£(z) 



z,--Bi(S(«)) 



£(z) 



y 6i (z)-r Ci (z) 

E(z) 



z 4 (1 - Bi(E(js)) 
Zi -Bi(E(z)) 



^(z)-^(z) Si(E(z))(l - Zi) 



Uz) 



(5) 



2i-B j (E(z)) 

A ? 



^ 1 - z« ^c,(z) (4) ^ Aj(l - Zj) cii 
2^ yJ7\ Z " ( 2^ E(z) cA h 



E(z) 7 t 

i=l v ; '* i=l 

which completes the proof. □ 
We would like to re-emphasise that =2 is completely determined when the are known. Thanks 
to the relations in Section 2.2, this is equivalent to knowing, for example, ^ . The latter functions 
were found for various polling models as mentioned in Section 1. 
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Remark 1. In Sidi et al. [10] the authors also focus on the joint queue length distribution at an 
arbitrary epoch. However, the result in [ ] is restricted to the case of non-zero switchover times 
and cycles with either all policies being gated or all being exhaustive. The authors averaged over 
N visit times and TV switchover times, but did not obtain our Theorem 1. Our theorem reveals that 
exactly the same structure holds, regardless of the service discipline. However, they do allow a 
more general customer behavior; their paper is one of the few polling studies in which the system 
is viewed as a network, with customers moving from queue to queue and the server visiting the 
queues cyclically. 

Remark 2. Observe that Theorem 1 immediately gives the formula for the marginal distributions. 
Indeed, for a vector zm,% = (1, ■ • ■ , 1, ^, 1, . . . , 1), £!(zm,i) = ^c{ z M,i)- From the 'step' (level 
crossing) argument it follows that S^ Ci [zu,i) is also the PGF of the queue length distribution in Qi 
at an arrival epoch at Qi. By PASTA it is also the steady-state distribution of Qi. 
Next take zt = (z, . . . , z). Theorem 1 now states that the PGF of the distribution of the total 
queue length (in terms of z) equals J2i=i ^i^ct ( z t)/ Ej=i ■ This formula may be interpreted as 
follows. By PASTA, J2(zt) is also the PGF of the distribution of the total queue length at an arrival 
epoch. By a level crossing argument, it follows that this equals the PGF of the distribution of the 
total queue length just after a departure epoch. The result now follows from the observation that 
a fraction A 2 / YLj=x or the departure epochs refers to a departure from Q,. 

Remark 3. It would be very interesting to have an interpretation of (10). Below we present a short 
derivation, which essentially combines several of the formulas above, but which fails to give an 
immediate explanation: 

JV N 
i=l i=l 

d2yi3) s% t (z) - Zi^M i 

k Pi ^u Mi - *i) EC k Ef =1 A,(i - *,) 
. XtiM^M - zi-rM) l El^jz) - v bi ( z )) 

m Eli Xi(y b ,(z) - Zt ,y Ci (z)) + Eti^M-^jz)) 
zZf=i x j( l - z j) E. v i -M 1 

Eill - z i) 

4. Workload distribution at an arbitrary time 

The idea of the proof of Theorem 1 can be also used to determine the LST W of the joint work- 
load distribution at an arbitrary epoch. For future use, let u; := (uj±, . . . , wjv) an d B(u>) : = 
(Si(wi), . . . , B N (u N )). Moreover, recall that £(z) = Y$=i A j(! - z j)- 

Theorem 2. For a general polling system as introduced in Section 2, let W(-) be the Laplace-Stieltjes 
transform of the joint workload distribution at an arbitrary time in steady-state. Then, 



(£(«)) -r^ (£(«)) 



lC^ x E(fl(u>)) E(B(w))- Wi 

where EC = s/(l — p)in the non-zero switchover case and EC = % 1 (0)/(A(l — p)) in the zero switchover 
case. 
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Proof The proof follows the same reasoning as the proof of Theorem 1. In particular, we focus 
only on the non-zero switchover times case. 

Let Xi(-) and Yi(-) be the LST of the joint workload distribution at an arbitrary moment during a 
visit to Qi and during a switchover (idle) period between Qi and Qi+i, respectively. Then by the 
stochastic mean value theorem 



N 

1 



Firstly note that using (13), 
Yi(u)=Yi(B(u)) 



1 ^(B(u,))-r 6i+1 (B(u;)) 



'< E(B(w)) 
Secondly, 

Xl (u) = x /°° /°° exp f-E(fl(«))« - w,t) dP (Sf ast < u, Br < t) . 

Bi(oJi) J u =oJt=o v ' 

Indeed, consider an arbitrary service time during a visit to Qi. The term S% i ^B{uy)/Bi{wi) cor- 
responds to the LST of the workload of all the customers in the system at the beginning of that 
service time, excluding the customer whose service is about to begin. Secondly, consider the work- 
load that arrives at all queues during the past part Bf as of that service time, and the residual part 
,B[ es of that same service time. Next integrate with respect to the joint distribution of Bf ast and 
BV es . Using (5) and the fact that, cf. [6, Section 1.6.3], 



f f exp (-E(B(w))u - Uit) dP (Bf ast < u, B[ es < i) 

Ju=0 Jt=0 ^ ' 



bt (e(B(o>)) - w, 



we obtain 

y bi (B(uj)) B^) - Bi(£(B(u))) li % i {B{u)) -% i {B{u)) B^i) - B^Biu))) 



Bi{ui) bihl{B{u))-Wi) b i Bi(ui) - Bi(E(B(w))) E(B(u>)) 
7i %,(B(uj)) — y Ct (B(u))) 



h E(B(u>)) - Wj 

This completes the proof. □ 

Remark 4. 5^ was found in [ ] for the more general class of polling systems with Levy input, but 
restricted to the case of branching-type service disciplines at all queues and non-zero switchover 
times. Notice that Y.(B(u>)) — uji is the Laplace exponent <f>f-(u)) of the Levy process 

Ai(t) = (Wi(t), . . . , Wi_i(t), - t, W i+ i(t), . . . , W N (t)), 

where Wi is a compound Poisson process with jump distribution Bi and rate A;. Moreover 
E(B(u;)) is the Laplace exponent of the Levy process W(t) = (Wi(t), . . . , Wjv(i)). After iden- 
tifying r 6i (B(a;)) and f Ci (B(uj)) of Theorem 2 with B^w) and B;(u>) of Theorem 3 of [5], it can 
now be verified that Theorem 2 indeed coincides with [ , Theorem 3]. 

Remark 5. Observe that Theorem 2 gives the LST of the distribution of W, the total workload in a 
polling system. With « T = (u,..., lu), we have E[e"" w ] = W{u) T ), so that 

- ^ - V ^(^))-^,(B(a; T )) 

s w-Efc=iAj(l-Bj(w)) ~[ Eli A 3 -(l- 4.(a;)) 
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Note that the LST of the amount of work Wm/g/i m a corresponding M/G/l queue, that is, an 



M/G/l queue with arrival rate A and service time LST E J= i XjBj(uj)/ A, is given by 

E[e~" Wjf/G/1 l = 



Note that Wm/g/\ is a ^ so m e amount of work in the polling system with zero switchover times. 
Finally, from the proof of Theorem 2 we know that the LST of the distribution of the total amount 
of work W aW itck in a polling system at an arbitrary epoch in a switching period is given by 

* (B(« r ))-^(B(« T )) 



(14) E[e~ uW '^] = -J2 sMut) = - 



Hence, we retrieve the well known decomposition (see, e.g., [3]) 

E[e~ uW ] = E[ e - uH ' 1 '/ G / I ]E[e- ulf » i "]. 

The distribution of W SW itch was only known in a few special cases, cf. [15]. From Theorem 2 it 
follows that in general its LST is given by (14). 
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